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We examine what happens when a strongly many body localized system is coupled to a weak 
heat bath, with both system and bath containing similar numbers of degrees of freedom. Previous 
investigations of localized systems coupled to baths operated in regimes where the back action of the 
system on the bath is negligible, and concluded that the bath generically thermalizes the system. In 
this work we show that when the system is strongly localized and the bath is only weakly ergodic, 
the system can instead localize the bath. We demonstrate this both in the limit of weak coupling 
between system and bath, and in the limit of strong coupling, and for two different types of ‘weak’ 
bath - baths which are close to an atomic limit, and baths which are close to a non-interacting 
limit. The existence of this ‘many body localization proximity effect’ indicates that many body 
localization is more robust than previously appreciated, and can not only survive coupling to a 
(weak) heat bath, but can even destroy the bath. 


Quantum localized systems violate many of the foun¬ 
dational assumptions of quantum statistical physics 
(such as the ergodic hypothesis), and present an exciting 
new frontier for research [T] . While localization was long 
believed to occur mainly in systems of non-interacting 
particles, the recent discovery of many body localiza¬ 
tion (MBL) [2]-[5] has ignited a blaze of interest in this 
field. It has been realized that quantum localized systems 
can display a cornucopia of exotic properties, including 
an emergent integrability IMS!, exotic quantum states 
of matter iMm, and unexpected behavior in linear [19] 
and non-linear m response. These properties not only 
dramatically revise our understanding of quantum statis¬ 
tical physics, but also offer a new route to dissipationless 
quantum technologies. A summary of progress in this 
field can be found in the review article [21] . 

Most works on MBL have focused on perfectly isolated 
quantum systems. Experimental systems, however, are 
always coupled (however weakly) to a thermalizing envi¬ 
ronment. The behavior of many body localized systems 
coupled to a thermalizing environment was first exam¬ 
ined in [221424] , in the limit where back action on the bath 
was negligible and the bath could be treated as Marko¬ 
vian. In this limit, it was argued that an arbitrarily weak 
coupling to a thermodynamically large bath should re¬ 
store ergodicity, thermalizing the system. These results 
suggested that perfect MBL would be unobservable in 
experiments, which would see instead only signatures of 
proximity to a localized phase. However, these works 
left open the question of whether different physics could 
result if the localization in the system were strong, and 
the bath were weak. Could many body localization then 
survive even after coupling to a heat bath? 

In this Letter, we show that when the system of interest 
is strongly localized, and the ‘heat bath’ is only weakly 
ergodic, then MBL in the system can not only survive 
coupling to the bath, but can even localize the heat bath. 
We call this phenomenon a ‘many body localization prox¬ 
imity effect,’ and it establishes that MBL is much more 
robust to coupling to an environment than was previ¬ 


ously appreciated. It also suggests a possible explanation 
for the numerical results recently presented in [261127], 
which counter-intuitively observed many body localiza¬ 
tion in an interacting model, when the non-interacting 
limit contained a single particle mobility edge. 

The system we consider consists of a D dimensional 
lattice which hosts two species of spinless fermions - c 
and d. The c fermions are present with density ric and 
have Hamiltonian 

iJc = E + UclciC^j + ^ Sicjci (1) 

(*J> » 

where tc is the hopping, U is a, nearest neighbor inter¬ 
action, and 5 is a random potential, drawn from a dis¬ 
tribution of width W. The width of the distribution is 
sufficiently large that the c particles in isolation are in 
an MBL phase, with a localization length We do as¬ 
sume that W is the largest scale in this Hamiltonian, and 
sets the characteristic energy scale of the MBL system. 
Meanwhile, the d particles are present with density rid^ 
and have Hamiltonian 

Hd = ^ tdd\dj + Xd\did\dj ( 2 ) 

{ij) 

This Hamiltonian in isolation describes a system in an 
ergodic phase (the ground state will be a Fermi liquid). 
The coupling between c and d systems is taken to have 
the form 

Hint ~ ^ ^ ( 3 ) 

i 

For simplicity in this work we assume ric = rid ^ although 
varying the ratio of the two densities would also be an in¬ 
teresting parameter to tune in future work. We consider 
two models. First, we consider a model (model I) where 
the coefficients gt are taken from some distribution (e.g. 
box distribution) of width G. Next, we consider a model 
(model H) where the gt are uniform, with strength G. 
Many of the results are easier to establish with a ran¬ 
dom coupling (model I), and we thus discuss this model 
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first, but we will subsequently show that a many body 
localization proximity effect can also arise with uniform 
couplings (model II). 

Weak coupling (small G): We first discuss the behavior 
with a weak random coupling gi taken from a distribu¬ 
tion of width G <^W. In this regime the disorder in the 
c sector is the primary source of randomness. A many 
body localization proximity effect can be shown to arise 
in two limits: the atomic limit ^ 0, and the non¬ 
interacting d limit A ^ 0. We establish by means of per¬ 
turbative expansions that this ‘many body localization 
proximity effect’ survives for small (but non-zero) and 
A respectively, such that there exists a finite sized region 
of parameter space where coupling a localized system to a 
bath results in localization of the bath. We subsequently 
show that this ‘many body localization proximity effect’ 
also arises if the coupling is uniform rather than random. 

Bath close to an atomic limit: t^ < G < W. We start 
by taking the limit td ^ 0. In this atomic limit for the d 
fermions, the eigenstates of the combined c and d system 
take the form 

i^)=i</>(w))®n4io) (4) 

fd 

i.e. d fermions are present on a certain set of sites {i}, and 
the c fermions are in a localized state, the precise wave 
function of which depends on the set of occupied sites {i} 
in the d sector. The wave function in the c sector depends 
on {i} because the distribution of d fermions affects the 
random potential seen by the c fermions. We note that 
a random distribution of atomic-limit d fermions will in¬ 
crease the effective random potential disorder in the c 
sector, from W to \/W‘^ -h and thus if the c fermions 
are localized at G = 0 (uncoupled c and d sectors), they 
should also be localized at non-zero G. 

We now turn on a non-zero hopping in the d fermion 
sector td 0^ and ask whether the resulting system is 
still localized. The hopping of a single d fermion changes 
the interaction energy between c and d fermion sectors by 
an amount of order G (we are working here with random 
couplings gi). It may also change the interaction energy 
in the d fermion sector by an amount of order A. Thus, 
naively a single hop by a d fermion takes the system ‘off 
shell’ by an random number drawn from a distribution 
of width at least G. However, the hopping of a single d 
fermion changes the potential seen by the c fermions, and 
thus in the c fermion sector acts like a ‘quantum quench.’ 
We now discuss to what extent the c fermion system may 
be able to relax to accommodate this change in energy. 

The c sector is in a localized phase, thus it will not be 
able to bring the system back precisely on shell. Rather, 
the c fermion sector should act like a ‘finite sized bath’, 
with the localization length setting the size, and s{T) 
being the entropy density. The relevant level spacing will 
thus be ^ IT exp(—s(T)^^). The hopping of a d fermion, 
combined with an appropriate relaxation in the c sector 


on length scales short compared to must thus take the 
system off shell by at least AE ^ W exp(—s(T)^^). 

To evaluate the convergence (or not) of the locator 
expansion [1], the change in energy must be compared 
to the associated matrix element. The matrix element 
between the two states is 

{{0\di 0 td4dj ((pij)) 0 d]\0)) « 

(5) 

Where and |0(j)) are the initial and final states in 
the c sector. We assume that the two states differ only 
on length scales less than (ignoring rare long range 
resonances). We make the standard assumption that 
the eigenstates are ‘ergodic’ on length scales less than 
i-e. that the initial state has similar overlap with all 
exp(s(T)^^) possible final states that differ from the ini¬ 
tial state only on length scales shorter than and also 
that the exp(5(T)^^) possible overlap matrix elements 
have random phases. Demanding normalization of the 
final wave function, we then conclude that a typical over¬ 
lap matrix element (0(i)|(/)(jf)) exp( —^^(T)^^) 

The locator expansion will converge if the typical hop¬ 
ping matrix element is smaller than the amount by which 
the hop takes the system off shell. From the estimates 
above, we conclude that a locator expansion in small td 
will converge iff 

td < inin(G, W exp(-is(T)^f)) (6) 

As long as this condition is satisfied, the combined c and 
d systems will both be localized, with eigenstates that are 
‘close’ to the form Thus in this there regime is a clear 
MBL proximity effect where coupling a localized system 
to a bath ends up localizing the bath instead. When this 
condition is violated, the locator expansion breaks down, 
which may indicate delocalization. 

The above analysis may be readily generalized to uni¬ 
form couplings of strength G (model II). As long as the 
density pattern in the c sector is spatially inhomogenous 
(which should be the case in the localized regime), the 
effective coupling between c and d sectors will be ran¬ 
dom, drawn from a distribution of width Gduc, where 
Suc is the width of the density distribution in the c sec¬ 
tor. An analogous argument may then be constructed for 
a many body localization proximity effect at weak td- We 
caution however that in the limit of site localization of c 
fermions tcfW 0, with uniform couplings, the effective 
disorder introduced in the d fermion sector is binary, and 
has percolating equipotential surfaces which break the 
locator expansion. This problem is absent with model I 
couplings. 

A weakly interacting low dimensional bath: A <C G ^ 
IT, td and D = 1, 2 

In this regime traditional locator expansions [1] fail, 
but we can still use weak localization results. Let us 
begin by taking A -> 0. In this limit, the d fermions are 
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described by a quadratic Hamiltonian, which takes the 
form 

Ha = Y,idd\dj + Y,Vid\di (7) 

{ij) * 

where Vi = gic\ci is a static random variable (since the c 
fermions are in a localized phase), and can be interpreted 
as a disorder potential, with the precise realization of the 
disorder depending on the state in which the c fermions 
are prepared. Thus, the d-fermion Hamiltonian describes 
fermions hopping in a random potential. It does not par¬ 
ticularly matter whether we work with model I couplings 
(random coupling constants) or model H (uniform cou¬ 
pling constants), since as long as the density pattern in 
the c sector is inhomogenous the d fermions will see an 
effective random potential anyway. We choose to work 
with model I for convenience and specificity. As is well 
known, free fermions hopping in a random potential will 
inevitably localize in D = 1, 2 mm, with a localization 
length that is power law large in td/G for D = 1 and 
exponentially large in td/G for D = 2. 

The argument above is cleanest in the site localiza¬ 
tion limit ^ O 5 when the c — d coupling is diag¬ 
onal in the c eigenbasis and there is no appreciable 
back action on the c system. At non-zero there 
will be an effective four d fermion interaction mediated 
by the c fermions. This four fermion interaction may 
be estimated in the manner of [16] and will be of or¬ 
der exp(— 2 /^c) for Cc < 1 and of order 

for > 1 (see supplement for de¬ 
tails). We are working here in units where the lattice 
scale has been set to one. Since an interaction is already 
induced by the c sector, we can also turn on the intrinsic 
interaction A 7 ^ 0. The matrix elements of this interac¬ 
tion in the basis of localized wave functions \ip) will be 
Aa/37<5 = ~ , ^here 

we have made use of normalization of the wavefunctions. 

The locator expansion should converge (such that the 
dressed eigenstates of the interacting problem are close 
to the eigenstates of the non-interacting problem), if the 
matrix elements of the interaction (both intrinsic and 
induced) are less than the accessible level spacing 
i.e. 

< 1; Cc < 1 

\wuC^ td) 

We note in particular that in one (two) dimensions, ^d 
is power-law (exponentially) large in td/G. For there to 
be a well controlled weak localization regime, we there¬ 
fore require not only that A is small, but also that either 
^ 1 or that W ^ td. We note also that the con¬ 
vergence criterion (i is simply an estimate obtained by 


considering the leading order terms in perturbation the¬ 
ory. A calculation to all orders along the lines of [3] is 
beyond the scope of the present work. 

As long as the interactions (both intrinsic and induced) 
are sufficiently weak, the wavefunction will take the form 

|$> = |0({a}))®n4(0)|O) (9) 

{«} 

i.e. the d fermions are in a product state where a set 
{a} of the non-interacting wave functions are occupied, 
the precise shape of the non-interacting wave functions 
depends on the state \(j)) in which the c fermions are pre¬ 
pared, and the c-fermions are in a state \(j)) which is many 
body localized, but depends in its detailed structure on 
the state of the d-fermions. 

Strong coupling (large G) We now point out that a 
many body localization proximity effect also arises for 
arbitrary A and t in the limit G ^ 00 . We demonstrate 
this for model H (uniform) couplings, but it is obviously 
true also for model I. 

In the large G limit, we have a description in terms of 
three species: ‘bound states’ whereby a c and d fermion 
sit on the same site, unbound cs, and unbound ds. A 
bound state cannot break apart or form, because this 
would change the energy by an amount of order G, which 
is the largest energy scale in the problem. 

The unbound c fermions are governed by a Hamilto¬ 
nian of the form Q, but now with certain lattice sites 
forbidden (the sites occupied by the c — d pairs and the 
unpaired ds). Since the c fermions localized with all lat¬ 
tice sites allowed, and the forbidding of certain lattice 
sites only presents an obstruction to transport, the un¬ 
paired c fermions should still be localized. Meanwhile, 
the bound c — d pairs are extremely heavy, with an effec¬ 
tive hopping matrix element of order tctd/G tc, and 
they see the same disorder potential as the c fermions. 
Thus, if the c fermions localized in a disorder potential 
of magnitude W, then the c — d composites should lo¬ 
calize also. Finally, that leaves the unpaired d fermions. 
These live on a random lattice, obtained from the origi¬ 
nal lattice by deleting all sites on which unpaired cs and 
c — d pairs are present. (This is because hopping onto 
an ‘occupied’ site changes the energy by an amount of 
order G, and in the large G limit this effectively makes 
occupied sites inaccessible). In one dimension, any finite 
density of deleted sites causes the ‘lattice’ on which the 
unpaired d fermions move to break up into finite sized 
segments, on which the d fermions are localized. Thus, 
in one dimension, in the G ^ 00 limit, we clearly have a 
system of three distinct simultaneously localized species, 
regardless of the values of td and A. 

In higher dimensions too, when a sufficiently large 
number of sites are deleted (large G limit), the lattice 
on which the unpaired d fermions live will break into dis¬ 
connected clusters, and at this point the d fermions will 
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most certainly be localized. Thus, there must be a regime 
in which the coupling of system and bath leads to the for¬ 
mation of localized species: localized c fermions, localized 
c — d composites, and d fermions which live on isolated 
islands that are separated by regions containing either 
c fermions or c — d composites. In fact, work on quan¬ 
tum percolation m suggests that for non-interacting d 
fermions (or, presumably, weakly interacting d-fermions), 
the localization transition precedes the percolation tran¬ 
sition, and d fermions will quantum localize even before 
the lattice on which they live breaks up into disconnected 
clusters. The density of deleted sites may be tuned by 
varying either the density of particles or the energy of 
the state (and hence the number of c — d pairs), and the 
quantum percolation (and hence localization) transition 
may also be tuned in this way. 

We note that in principle there is an ‘exchange’ process 
by which an unbound c and d^c — d composite on adjacent 
sites could exchange position by moving the d fermion, 
without taking the system off shell. This is really only 
possible if the c fermions are in the site-localization limit 
W ^ tc (since the c — d composites are always site local¬ 
ized in the large G limit). However, in this limit of site 
localization for c fermions, localization requires the addi¬ 
tional condition that the set of possible energy conserving 
d exchanges should not percolate. 

If G is large but not infinite, then d fermions will in 
principle be able to ‘tunnel through’ the barriers pre¬ 
sented by c fermions, but virtual processes involving c — d 
composites will then also give rise to a random onsite 
potential for the lattice on which the d fermions live (in¬ 
herited from the inhomogenous density pattern for the c 
fermions). In low dimensions and in the non-interacting 
limit A —> 0, and for sufficiently large G, this will be suffi¬ 
cient to localize the d fermions, and arguments along the 
lines of [3] suggest the d fermions will continue to be lo¬ 
calized for small but non-zero A. Of course, the existence 
of a ‘many body localization proximity effect in the limit 
of weakly interacting d fermions is not particularly sur¬ 
prising, since something similar occurs for weak G. The 
truly novel feature is the existence of an MBL proximity 
effect at arbitrary A in the limit G ^ oo. 

Conclusion: Thus, we have demonstrated that when 
a many body localized system is coupled to a bath, the 
end result can be the survival of MBL in the system, and 
the localization of the bath. We have demonstrated that 
this happens in three distinct limits: a bath close to the 
atomic limit that is weakly coupled to a strongly local¬ 
ized system, a weakly interacting bath weakly coupled to 
a strongly localized system, and an bath (with arbitrary 
bath parameters) that is strongly coupled to a localized 
system. The problem of localized systems coupled to 
baths is thus much richer when the back action on the 
bath is taken into account, and indeed the coupling of a 
localized system and a delocalized system can result in 
the localization of both. This may explain recent numer¬ 


ical results PH HI]. The development of a general theory 
of localized systems coupled to baths is an interesting 
challenge for future work. 
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order in perturbation theory in weak G. Solid lines denote d fermions and dashed lines denote c fermions. The state of the d 
fermions is labelled by Greek characters, while the state of the c fermions is labelled by English characters. We assume the 
localization centers of the d fermion wavefunctions (l)cx,/ 3 ,'y,s are all within of each other. 
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SUPPLEMENT: EFFECTIVE FOUR d FERMION INTERACTION MEDIATED BY c FERMIONS 


When ^ 0, the coupling of c and d fermion sectors leads to an effective four d fermion interaction. In this 
supplement we estimate this interaction, first for < 1 and then for > 1- The calculation follows [16]. The two 
processes contributing to the effective interaction at lowest order in perturbation theory in weak G are shown in Fig{^ 
The process (i) dominates and will be discussed here, but the contribution of process (ii) can be calculated similarly. 
The process Figj^i) represents an effective T matrix element of the form 


pabc 


i 


= G" / d^nd^r2 




C( ri)0/3(ri)07(r2)</>|(r2)^a(ri)V>b(ri)^^(r2)V>c(r2) 
AE 


( 10 ) 


where (p represents a c fermion wavefunction, pj represents a d fermion wavefunction, and AE represents the amount 
by which the intermediate state is off shell. 

If < 1 then the above matrix element falls off exponentially in |ri — r 2 | and the integral is effectively restricted 
to neighboring sites. Meanwhile, the scattering in the c fermion sector takes us off shell by an amount of order 
W. Finally, since each site has in effect a single associated state 'ipa,b,c^ the matrix element is of order « 

(G^/VF)^^^^ exp(— 2 /^c) exp(i 0 abc )5 where we have assumed the d wavefunctions p are localized to a volume of radius 
and where pabc is a random phase. The effective induced interaction may be obtained by summing over all 
intermediate c fermion states i.e. 


\^f f \ ^ rjiabc 

abc 


( 11 ) 


There are terms contributing to the above sum (since the four d fermions involved can interact through c fermions 
anywhere in the volume where they all overlap). Each term has similar magnitude but random phase. Summing over 
all the terms in rms fashion gives an additional factor of ^ , and as a result we conclude that the effective four 
fermion interaction mediated by the c fermions is of order (G^/VF)^^ ^d /2 as quoted in the main text. 

The calculation is similar when > 1- The main differences are (1) the integrals / d^rid^r 2 now go over volumes 
of linear dimension (2) normalization requires that \%Ij\ ^ within of the corresponding localization center 

(3) the minimum energy change is AE ^ and (4) instead of having terms contributing to the sum (11), 

there are contributing terms (i.e. the volume can be divided into (Cd/Cc)^ cells, in each of which there are 

possible choices of a, 5, c). As a result, the effective four fermion interaction becomes of order (G^/VF)^^ 
again as quoted in the main text. 







